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ABSTRACT 

The normalized non-dimensional von Karman-Howarth 
equation for isotropic homogeneous decaying and forced 
steady turbulence is integrated to obtain expressions for the 
dissipation rate coefficient C £ = (Le)/(i< 2 ) 3 / 2 , where L de- 
notes the longitudinal integral length scale, e the mean dis- 
sipation rate and (u 2 ) the mean variance of the longitudinal 
velocity fluctuations. For decaying turbulence the final exact 
expressions for C £ for the low and high Reynolds number limit 
depend on the decay exponent n, which is known to depend 
on the initial velocity structure at the turbulence production. 
The dependence on n leads to a non-universal coefficient. The 
expressions for the steady forced case depend on the forcing 
mechanism and thus are not universal either. Nonetheless, a 
lower value and considerably less scatter as compared to the 
decaying turbulence case should be expected when similiar 
forcing algorithms are employed. 



INTRODUCTION 

A fundamental assumption on which the phenomenol- 
ogy of turbulence is based, is the independence of the mean 
energy dissipation e of the viscosity v, provided the latter is 
very small. This concept, also termed dissipation anomaly, cf. 
Frisch 1 1 1, has first been introduced by Taylor |2| in 1935 and 
it directly translates, based on dimensional grounds, to a pro- 
portionality between e and (m 2 ) 3 ' 2 /L with L being an integral 
length scale. The assumption of a universal proportionality 
constant plays a fundamental role in turbulence modeling and 
even directly determines one of the model constants in today's 
widely used one- and two-equation RANS models, such as the 
k — e model or the k — (0 model, where the unknown Reynolds 
stresses are closed based on an eddy viscosity ansatz. How- 
ever, as early as 1953 when Batchelor |3| plotted values of 
the proportionality constant in his textbook the question arose, 
whether or not it is truly universal and if it is not, how it de- 
pends on the flow configuration under consideration. Exper- 
imental data from various types of decaying turbulent flows 
such as wakes, jets and grid turbulence (with different geome- 
tries) have been analyzed and the resulting values of C £ have 
been collected for instance by Sreenivasan (4]|5] and Burat- 
tini (6). The large scatter of the data, even for similar flow 



configurations, has led all of the above mentioned authors to 
the conclusion that the constant is not universal and depends 
critically on the way the turbulence is generated. Recently 
Goto and Vassilicos [7 J showed that the value of C £ depends 
on the internal stagnation point structure and can thus not be 
universal, as the stagnation points themselves depend on the 
large scales of the flow. The purpose of this work is to de- 
rive expressions for C £ based on the von Karman-Howarth 
equation, for the case of homogeneous isotropic decaying tur- 
bulence and homogeneous isotropic forced steady turbulence. 
The paper is organized as follows. In Chapter 1 we will de- 
rive exact expressions for the dissipation rate coefficient for 
decaying and forced turbulence based on the normalized von 
Karman-Howarth equation. The results of the analysis are dis- 
cussed and compared with each other in chapter 2. In chapter 
3 a distinction is made between the limit of large and small 
Reynolds numbers and exact and approximate results for the 
dissipation rate coefficient are derived for both cases. Con- 
cluding remarks are given in chapter 4. 



1 ANALYSIS 

For isotropic homogeneous turbulence we define the 
non-dimensionalized two-point longitudinal second and third 
order structure functions as 
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where u(x,t) denotes the instantaneous longitudinal ve- 
locity fluctuation at point x which is seperated by a dis- 
tance r from the second point. The non-dimensional sec- 
ond order structure function can be transformed into the non- 
dimensional two-point correlation yielding 
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where angular brackets denote ensemble averages. 
Let us define two length scales 



L(t)= f°f{r,t)dr 
Jo 
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the ratio of which is precisely the above introduced en- 
ergy dissipation rate coefficient as 



(the superscript d stands for decaying) with Rei = 
( (u 2 ) l ' 2 L) I v a Reynolds number based on the integral length 
scale defined in eq. [3] 

Noting that by definition JJ° fdf = 1 we integrate eq. [5] 
to obtain 
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1.1 Decaying turbulence 

We write the von Karman-Howarth equation, cf. von 
Karman and Howarth [8], following an approach by Lund- 
gren |9J as 
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where the second integral on the left hand side is ob- 
tained by partial integration. Eq. fTOlfhus simplifies to 
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where v denotes the molecular viscosity of the fluid. 
Let us define the new coordinates 



f = r/L{t), T = t. (6) 
Then the time derivative df(f, t) /dt in eq. |5]becomes 
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The derivative df/dt may be expressed as 
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Using relation|4]we write L = C £ (u 2 )^l 2 / £ so that 
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with du 2 /dt = — 2/3e for isotropic decaying turbulence. 
Under the assumption that the turbulence has reached an 
asymptotic state at which the decay of the turbulent kinetic 
energy follows a power law decay of the form k(t) °c t~" 
with a constant decay exponent n, which corresponds to self- 
preserving solutions of the von Karman-Howarth equation 
(for a detailed discussion see Dry den |10| , Batchelor |l 1 1 . Ko- 
raeyev 1 12]), one obtains 
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where T; = k/e defines an integral eddy turn over time 
scale. Eq. 13 can be further simplified as T; = t/n and 
dlne/dt = - (n+l)/f, so that 



Introducing the new coordinates into eq. |5](and dropping 
the functions' arguments for simplicity) yields 
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Eq. [T5] constitutes an exact relation for homogeneous 
isotropic decaying turbulence. At this point it is pertinent 
to make reference to Rotta's [13] work, who used Heisen- 
berg's [ 14J theory in combination with Loitsianskii's 1 15| in- 
variant hypothesis to derive explicit expressions for for the 
low and high Reynolds number limit of decaying turbulence. 
Different from Rotta's approach, the above analysis does not 
assume any invariant integral of motion or ad-hoc assumptions 
and the final exact expression (eq. \15\ preserves the explicit 
dependence on the decay exponent. 



The integral on the right hand side of eq. [20] can be 
solved analytically 




We thus obtain the equivalent expression of eq. |15|for 
the forced case, 



1.2 Forced turbulence 

In the case of forced turbulence we write the * L ( f°° (—djix) 12 [°° _i df 
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von Karman-Howarth equation following an approach by E If \Jo ? Rei Jo df 

Fukayama 1 16) as 
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where P denotes the mean energy input through the forc- 
ing and the last term on the right hand side of eq. |T6] is the 
forcing- velocity correlation term with 
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Here j is a dummy variable and the forcing is assumed to 
be random Gaussian, delta correlated in time, homogeneous 
and isotropic. Different from Fukayama [ 16] we assume the 
forcing to have compact support limited to the largest scale of 
the system denoted by If. Eq. [T6]is thus valid for < r < If. 
For a detailed derivation see Fukayama 1 1 6| . In the steady 
case we have 



dD, 



dt 



= 0, /> = £, 



(18) 



and the mean energy input must equal the mean dissi- 
pation. We write eq. [T6]in non-dimensional form using the 
coordinates defined in relation|6] so that we obtain 
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Integration over the entire domain (where we implicitely 
let U ■ — > 00) yields 
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2 DISCUSSION 

Eq. [T3]and eq. [22] constitute two exact results for the 
dissipation rate coefficient in homogeneous isotropic decay- 
ing and forced steady turbulence, respectively. A direct com- 
parison of the two expressions reveals that they are identical 
in form and depend on the same integral expressions involv- 
ing the non-dimensional third order structure function drjj.(f) 
and the non-dimensional two-point correlation f(f). These 
integral expressions are only multiplied by a different factor, 
which in the case of decaying turbulence depends on the decay 
exponent n and in the case of forced turbulence on the ratio of 
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the longitudinal length scale L, defined in eq. [3]to the scale 
at which energy is pumped into the system to sustain a steady 
turbulent flow. For the following discussion we assume the 
high Reynolds number limit of eq. [15] and eq. [22] for which 
the last term on the right hand side proportional to Re^ 1 can 
be neglected. In this limit the dissipation rate coefficient will 
be denoted as C e ,c» and its expression only involves an integral 
over the non-dimensional third order structure function whose 
functional form is a-priori not known. As the upper integra- 
tion limit extends to infinity one can deduce that the integrals 
involve contribution from the largest scales of the turbulent 
field which cannot be expected to be universal. For the case 
of decaying turbulence the decay exponent dependent factor 
can be expected to be different for different types of flows, 
as the decay exponent n which is realized for the flow de- 



pends on the initial conditions as stated by George |17| 
and Barenblatt [19|. Experimental results from Ling 
have supported the non-universality of the decay exponent. 
This leads to the conclusion that the energy dissipation coef- 
ficient in the high Reynolds number limit cannot be universal 
for homogeneous isotropic decaying turbulence. In the case 
of forced turbulence the value of the factor L/lf cannot be 
specified from first principles. Nonetheless, one can draw the 
valuable conclusion that it must be smaller than the decay ex- 
ponent dependent term for decaying turbulence. Assuming 
If to be of the size of the system (the computational domain 
for DNS calculations) then, as the non-dimensional two-point 
correlation f(f) decays to zero for r—¥lf and /(0) = 1, eq. 
pleads to the conclusion that L < If, so that L/lf < 1, while 
for all possible values of the decay exponent, the coefficient 
In I (2>n — 2) > 1. One must thus conclude that in the high 
Reynolds number limit the energy dissipation coefficient as- 
sumes smaller values for forced turbulence as compared to 
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decaying turbulence. In addition the lack of a strongly vary- 
ing parameter such as the decay exponent n will lead to less 
scatter in the case of forced turbulence. 



3 EXACT AND APPROXIMATE SOLUTIONS 
3.1 Low Reynolds number limit 

For the low Reynolds number limit the dissipation rate 
coefficient is denoted by Cf . This regime corresponds to the 
viscous dominated final period of decay which is also char- 
acterized by a power law decay and during which the contri- 
bution of the third order structure function in eq. [5] can be 
neglected. The decay exponent however assumes values dif- 
ferent from the ones of the high Reynolds number regime. As 
the von Karman-Howarth equation is for this case in closed 
form, its general solution can be expressed in terms of a spe- 
cial function, the confluent hypergeometric function to be de- 
noted by M(a,P,z), where a and j8 are parameters and z is 
the coordinate. The two-point correlation can then be written 
as 
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where integration by parts and the fact that t = 
3/2n (u) 2 /e and Re l ^ 2 = C^J 2 Re L have been used. The 
second integral on the right hand side is more straightforward 
to manipulate and yields 

(29) 

The final expression for the low Reynolds number limit 
reads 



u 2 f = u 2 M(ii,5/2,-$ 2 /S), 



(23) 



'£,0 " 



3n ( 


rr 


[2n\ 




2-2,11 y 


V 6n 


VtJ 


L 



Rel 



(30) 



where n is the decay exponent and § = r/ y v(t — to) a 
similarity variable, to is a virtual time which can be neglect 
against t as the low Reynolds number limit in decaying turbu- 
lence corresponds to large times t. For a detailed discussion 
of the final period of decay see Barenblatt |19| . 

Differentiating eq. [23]wifh respect to / yields 
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Dividing eq. u4|by (« 2 ) 3 / 2 and integrating yields, 
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The left hand side of eq. [25]can be simplified to 
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Eq. [30]is an exact expression for the final period of de- 
cay in homogeneous isotropic decaying turbulence. It obeys 
the well-known scaling C £ q Re^ 1 while the proportionality 
constant depends on the decay exponent n during the final pe- 
riod of decay (the integral over the confluent hypergeometric 
function is also fully determined by the value of n). 

Eq. [30] constitutes an exact relation for which without 
further assumptions numerical results can be obtained and it 
reveals an explicite dependence on the decay exponent n that 
governs this regime. The value of n however depends, as it 
does in the high Reynolds number limit on the initial condi- 
tions of the turbulence production, cf. Barenblatt |19| . As 
a consequence, the behaviour in the low Reynolds number 
limit cannot be expected to be universal. As stated by Baren- 
blatt 1 1 9 1 , solutions to the problem with a value of n > 5 /2 are 
not of physical interest as they are structurally unstable with 
respect to the initial conditions while the range of possible 
values for 1 < n < 5 12 is continuous. The coefficient of the 
right hand side of eq. 30 in front of Re^ 1 can be calculated for 
a given value of n. Table[TJshows this value for different val- 
ues of n. Note that for n = 5/2, which corresponds to the so 
called Loitsianskii case, the confluent hypergeometric func- 
tion reduces to a Gaussian error function so that the integral 
over it can be expressed analytically as 



The first integral on the right hand side of eq. [25] can be 
simplified with 
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The values obtained here are of the same order of mag- 
nitude as the ones reported by Bos |21| who derived an ex- 
pression for the dissipation rate coefficient based on modelled 
energy spectra. For the final period of decay he finds for 
Cg ~ 19..32We^' depending on the assumptions made on 
the energy spectrum. 
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Table 1 . Proportionality coefficient C d £ () <=c Re L 1 for the 
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period of decay. 



The expression for the dissipation rate coefficient in the 
high Reynolds number limit for forced steady turbulence reads 
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3.2 High Reynolds number limit 

Different from the low Reynolds number limit, no nu- 
merical values for the energy dissipation rate coefficient can 
be deduced from eqs. [15] or [22] without further knowledge of 
the form of dux(f). One has thus to resort to a closure ap- 
proximation for diLi(f) to solve the corresponding integrals. 
We will employ an eddy viscosity closure formulated in real 
space as proposed by Oberlack and Peters |22| of the form 
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In eq. 
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the constant a = 6/(5C 3 ' 2 ) has been cho- 
sen to correctly reproduce Kolmorogov's four-fifth-law for the 
third order structure function Dm = 4/5 (er) and the two- 
third scaling for the second order structure function Dn = 

2/3 

C (er) ' for small values of r, where C is Kolmogorov's con- 
stant. Caution has to be exercised as such a closure assump- 
tion neglects the non-universal influence of the large scales on 
the functional form of djjjjj) and thus only one single value 
for the integral can be predict independent of the flow type. 

The integration of the right hand side of eq. [T5]and eq. 
[22]yields 
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The expression for the dissipation rate coefficient in the 
high Reynolds number limit for decaying turbulence then 
reads 
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Sreenivasan 1 23] gives a compilation of values of Kolo- 
mogorov's constant C taken from different experiments and 
finds C = 2.0 ±0.4. Choosing these limits and varying the 
decay exponent n in a range from 1.0 to 1.9 corresponding to 
the values observed in experiments and DNS calculations (cf. 
Ling and Wan |20| . Uberoi [24| and Corrsin [25] and refer- 
ences therein), one obtains values of C{? «, ranging from 0.6 to 
2.2, which fall in the same range as the values quoted in the 
literature (cf. Sreenivasan (4][5j and Burattini [6]) and repre- 
sent a large part of the scatter observed. 



Different from the theory in the DNS calculations energy 
is pumped in over a finite band of large scales and not at a 
single scale If. To account for this difference, an equivalent, 
single wavenumber )Cy is calculated, that corresponds to the 
band of forced scales < K < Kq as 



_ f*°KE{K)dK _ q+l 
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where the energy spectrum at the large scales is assumed 
to follow E(k) oc K q , cf. Goto jT). For a calculation domain 
of size 2n this yields 
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Note that the approach taken in eq. [36]is only one possi- 
ble way of determining a single forced scale and corresponds 
to calculating the center of gravity of the energy spectrum in 
the range If < r < °°. 

The ratio L/lf can for example be determined from the 
calculations carried out by Goto |7 | and it turns out that it 
varies slightly between 0.42 and 0.49, confirming the hypothe- 
sis that L<lf. With these values and varying the Kolmogorov 

constant in the above ranee the interval 0.28 < c£„ < 0.60 
can be deduced. 

Let us briefly discuss the above found values for C E ^ in 
the context of RANS turbulence modeling. One common ap- 
proach to close the unknown Reynolds stresses in RANS tur- 
bulence models is to employ an eddy viscosity closure (26), 
where assuming the energy dissipation coefficient to be uni- 
versal in the high Reynolds number limit, one can show that 
the model constant Cn is related to the energy dissipation co- 
efficient, yielding Q.oo = (3/2) 3/2 c^ 4 . 

With the commmon value of Cji = 0.09 this leads to a 
value of C e .oo = 0.30 which is close to the lower bound pre- 
dicted for forced turbulence and considerably lower than any 
value predicted for decaying turbulence. In the light of the 
non-universality of the energy dissipation coefficient the rela- 
tion of the model constant to the theoretical energy dissipation 
coefficient is however rather hypothetical and the model con- 
stant cannot be predicted from first principles but rather an 
optimal value for different flow configuration has to be cho- 
sen. 



4 CONCLUSION 

An analysis based on the integration of the von Karman- 
Howarth equation in the case of isotropic homogeneous de- 
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caying and steady forced turbulence has led to three exact re- 
lations for the energy dissipation rate coefficient C £ . Of spe- 
cial interest is the high Reynolds number limit of these ex- 
pressions and the question whether the energy dissipation rate 
coefficient attains a universal value for infinite Reynolds num- 
bers. It could be concluded, that in the case of decaying turbu- 
lence no universal constant is attained as the final expression 
retains an explicite dependence on the decay exponent of the 
flow which is known to be non-universal. As for the steady 
forced case it is found that the final expression for C £ depends 
on the forcing mechanism. Although different forcing mecha- 
nisms are employed by different authors, it is concluded that, 
the values obtained from such DNS calculations will be close 
to each other when similar forcing schemes are employed and 
that they exhibit less scatter than in decaying turbulence, as a 
strongly varying parameter is absent. Based on the exact rela- 
tions it could be concluded that the energy dissipation rate co- 
efficient is smaller for forced turbulence that for decaying tur- 
bulence. Numerical values however cannot be obtained from 
first principles only. The introduction of a closure assumption 
for the third-order structure function allows the calculation of 
numerical values for both cases. These turn out to depend on 
the value of the Kolmogorov constant. For both decaying and 
forced turbulence a corridor of numerical values can be given 
which agrees fairly well with the values reported in the litera- 
ture. 
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